Violation of local realism in a high-dimensional two-photon setup with non-integer 

spiral phase plates 
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We propose a novel setup to investigate the quantum non-locality of orbital angular momentum 
states living in a high-dimensional Hilbert space. We incorporate non-integer spiral phase plates in 
spatial analyzers, enabling us to use only two detectors. The resulting setup is somewhat reminis- 
cent of that used to measure polarization entanglement. However, the two- photon states that are 
produced, are not confined to a 2 x 2-dimensional Hilbert space, and the setup allows the probing of 
correlations in a high-dimensional space. For the special case of half-integer spiral phase plates, we 
predict a violation of the Clauser-Horne-Shimony-Holt version of the Bell inequality (S < 2), that 
is even stronger than achievable for two qubits (S = 2\/2), namely S = 3^. 

PACS numbers: 03.67.Mn, 42.50.Dv 



I. INTRODUCTION 

Recently, the orbital angular momentum (OAM) of 
light has drawn considerable interest in the context of 
quantum information processing. The spatial degrees of 
freedom involved in OAM 1] provide a high-dimensional 
alphabet to quantum information processing (i.e. quAits 
instead of qubits) @,0|. Additionally, since OAM is asso- 
ciated with the topology of the electromagnetic field, the 
use of this observable in quantum entanglement may lead 
to states that are inherently robust against decoherence 
1. 

The most popular OAM analyzer when dealing with 
conservation, correlation and entanglement of OAM con- 
sists of a so-called fork hologram Q, i.e. a binary phase 
hologram containing a fork in its center 6] , together with 
a spatial-mode detector consisting of a single-mode fiber 
connected to a single-photon detector; such analyzers 
have been used in the three-dimensional case, i.e. N = 3, 
by Vaziri et ai, 3]. In that experiment, proof of entan- 
glement of the OAM degree of freedom of two photons 
was given by showing that a generalized Bell inequality 
was violated; this scheme requires 6 detectors, namely 3 
in each arm, and one has to measure 3x3 coincidence 
count rates to perform a measurement for a single 
setting of the analyzers. 

In the present paper, we consider the use of spiral 
phase plates (SPPs) Q instead of phase holograms in 
an OAM entanglement setup, enabling us to investigate 
high-dimensional entanglement with only two detectors. 
More specifically, we will consider SPPs that impose on 
an optical beam a non-integer OAM expectation value 
per photon, in units of h Q- With such devices, com- 
bined with single-mode fibers to form quantum-state an- 
alyzers, we propose to build an OAM entanglement setup 
that is reminiscent of the usual setup to measure polar- 
ization entanglement 8] , where the rotational settings of 
the analyzers (polarizers in that case) is varied. We will 
show that it is possible to identify SPP analyzer settings, 
in the spirit of horizontal and vertical aligned polarizers, 



when using half-integer SPPs, allowing observation of 
high-dimensional entanglement (A > 2), in contrast to 
the polarization case (N = 2). These claims are sup- 
ported by calculations; we predict highly non-classical 
quantum correlations (S — 3^), i.e. stronger quantum 
correlation between two photons than the maximum cor- 
relation between two qubits (S = 2y/2). 



II. SPIRAL PHASE PLATES 

A SPP, shown in Fig. Qfa), is a transparent dielec- 
tric plate with a thickness that varies as a smooth ramp, 
thus phase shifting an incident field linearly with the az- 
imuthal angle 9 pj. As can be seen, the plate carries a 
screw discontinuity, expressed by the spiraling thickness, 
and an edge discontinuity, i.e. a radially oriented step 
with height h s . The difference between the maximum 
and minimum phase shift is written as 27r^, where t is not 
necessarily integer; in fact, t depends on depends on the 
step height h s , the difference in refractive indices of the 
SPP and surrounding medium, and the wavelength of the 
incident light 0. Thus, a photon propagating through 
this plate will acquire an OAM with expectation value 
equal to £H 00. Placing such a plate in the waist of 
Laguerre-Gaussian beam, the field in the exit plane just 
behind the plate will be described by 



(r,0\S(e)\l,p)=u l l G (r,0)exp(ie0) 



(1) 



where \l,p) are the Laguerre-Gaussian (LG) field states 
and S(£) is the operator representing the effect of the 
SPP on the input mode. Note that we have neglected 
the uniform phase shift that is caused by the SPP's base 
with height ho , since it acts as a plane-parallel plate. By 
adding another plane-parallel plate with the appropriate 
thickness, the total phase shift of these two can be made 
equal to an integer multiple of 2ir. 

The function u l ^ G (r, 0) in Eq. is the complex ampli- 
tude of the Laguerre-Gaussian beam in its waist plane, 
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FIG. 1: (a) Schematic drawing of a SPP. The device shifts 
the phase of an incident beam proportional to the azimuthal 
angle 0. (b) A calculated far-field diffraction pattern of a 
fundamental Gaussian beam after propagating through an £ = 
3^ plate positioned in its waist plane, showing that rotational 
symmetry is broken. Black and white denote low and high 
intensity, respectively. 

given by 

(r,6\l,p) =uf p G (r,6) 

( r2 \ 
x exp exp (U9) , 

V O 

where w§ is the waist radius, L l p (x) an associated La- 
guerre polynomial |10| and Ci p a normalization constant. 
In the paraxial limit, the LG free-space modes, enumer- 
ated by the (integer) indices I and p, form a complete 
basis of spatial modes. The LG index I, which should 
not be confused with the SPP step index £, is related to 
the OAM that is carried in that LG mode, namely lh per 
photon 0,013, while the index p provides information 
on the number of nodes along a transverse radius of the 
mode. When placing a SPP in a beam in a single LG 
mode, the output will generally be in a superposition of 
LG modes and will be no longer invariant under free- 
space propagation. For integer values of the step index 
£ of the SPP, the edge discontinuity is effectively absent 
and this super pos ition will only be with respect to the 
LG-index p fA lllf. In that case, the intensity distribu- 
tion of the mode will be doughnut-shaped in the far field. 
For non-integer £ values, the edge discontinuity does not 
vanish and the superposition of modes will also be with 
respect to the index /, which is related to the OAM of 
the mode 0. Such superpositions consist in principle of 
an infinite number of LG components. Effectively, this 
number is finite and increases with l\ as an example, if 
£ = \, 11 LG components are sufficient to describe 87% 
of the field behind the SPP, while for I = §, 224 LG 
components are required. As SPPs with non-integer £ 
can create such high-dimensional superpositions of OAM 
modes, we anticipate that, when employed suitably, they 
can also project onto such superpositions. We therefore 
propose to incorporate such non-integer SPPs in analyz- 
ers for OAM states living in a high-dimensional Hilbcrt 
space. 

From a topological point of view, SPPs with non- 
integer £ imprint a mixed screw-edge dislocation on an 



incident field. The result is rotational asymmetry of the 
imprinted phase distribution and thus of the emerging 
field, which becomes visible in the far-field intensity pro- 
file (Fig. Hfb)). K is the orientation of the step in the 
transverse plane, that we wish to exploit as an analyzer 
setting in a new bipartite entanglement scheme. 

Since £ shall be chosen to have a non-integer value, it 
is important to realize that, when an incident field passes 
through an SPP in combination with its complement (i.e. 
a SPP with the same step height and orientation, but an 
inverted vorticity), the beam basically passes through a 
plane-parallel plate that shifts the phase of the field, in 
an azimuthally uniform way, by 2ir£, 

^compi (i)S{l) = exp (i 2tt£) I, (3) 

where / is the identity operator and where we keep the 
exponent since exp (i2n£) ^ 1 for non-integer I. Since 
S(£) is unitary, it follows that 

^compl ^ = exp £ t (£) f (4) 

where S'(l) is the Hermitean conjugate of S(£). As S^(£) 
and S compl (£) only differ by a multiplicative phase fac- 
tor we can again use a plane-parallel plate in the exper- 
iment to compensate for this phase factor. Similar to 
the uniform phase shift caused by the SPP's base ho, we 
will neglect this phase shift as well, as it can be trivially 
dealt with. Henceforth, the operator for the compensat- 
ing SPP, with inverted vorticity, will be represented by 
the operator S' {£) . 

III. PROPOSAL FOR AN EXPERIMENT 

In the experiment on OAM qutrits (N = 3), fork holo- 
grams were used @- Those holograms can only modify 
the OAM expectation value by an integer number |l2j, 
depending on the diffraction order of the hologram. It 
also required the use of three analyzers in each arm of a 
spontaneous parametric down-conversion (SPDC) setup. 
In contrast, our proposed OAM entanglement experiment 
uses non-integer SPP state detectors and only requires a 
single analyzer-detector combination in each arm. This 
scheme is shown in Fig. El it has been inspired by the 
setup used in polarization entanglement The SPPs 
that are inserted in the two arms should be chosen to 
obey conservation of OAM (see also section IVT1 . Thus, 
if the classical pump beam does not contain any OAM, 
a typical choice for the signal and idler SPPs would be 
to have step indices I and —£ respectively, where, in our 
case, £ should have a non-integer value. The role of the 
lenses in signal and idler path, L s and Lj respectively, will 
be discussed in section lVll Finally, single-mode fibers are 
indicated by F s and F;. 

By manipulating the transverse axes of our analyzers, 
namely the SPP steps, we have access to various pho- 
ton states that live in a high-dimensional OAM Hilbert 
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FIG. 2: Proposed experimental setup. A nonlinear crystal 
(NLC) splits a pump photon in a signal photon and an idler 
photon by the process of SPDC. In each path, a SPP (SPP S , ; ) 
is inserted with a single-mode fiber F B ,i, together forming the 
analyzer. The coincidence count rate of detectors D s ,i is mea- 
sured as a function of the SPP angular settings a s and «i. 



space, as mentioned earlier. We may orient these edges 
arbitrarily in the transverse plane, as shown in Fig. 
thus allowing their use as angular analyzers. When com- 
bined with a single-mode fiber, the SPP with step index 
£, set at an azimuthal angle a s , projects the incident 
photon state out onto the OAM state with expectation 
value — I with edge angle a s . As we will argue in sec- 
tion IVII the coincidence count rate will depend, like in 
polarization entanglement, only on the relative angle of 
the transverse axes. 

We stress that, in spite of the superficial similarity 
between a polarizer and a non-integer SPP, they are of 
course very different devices; for example, while polariza- 
tion corresponds to alignment, the SPP edge corresponds 
to orientation. In other words, with a polarizer, one can 
analyze the alignment of the electrical field oscillation, 
either horizontal or vertical, thus yielding a periodicity 
of 7r when rotated. With a SPP with non-integer £ one 
can analyze the spatial orientation of a field, thus yield- 
ing a periodicity of 2ir when rotated (see e.g. Fig.QIb)). 
Thus we expect that the coincidence count rate will have 
a periodicity of 2n when one of our analyzers is rotated. 

An equally important difference between the two cases 
is that, whereas polarization Hilbert space is two dimen- 
sional, OAM Hilbert space is infinite dimensional. As 
we will see, this makes rotation of SPPs fundamentally 
different from rotation of a polarizer. In order to ad- 
dress these aspects explicitly, we need a basis of the OAM 
Hilbert space that is suited for our purpose. 



IV. NON-INTEGER OAM STATES 

Our aim is to construct a complete, orthonormal basis 
that contains non-integer OAM states as basis elements. 
To this end, we consider all states in the polar represen- 
tation and separate the radial and angular parts, so that 
the arbitrary state \ f) = \r,9) can be written as |r)|r), 
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Azimuthal angle 9 

FIG. 3: Plot of the phase shift imprinted by a SPP as a 
function of the azimuthal angle. The solid line shows this 
shift for a plate with its edge oriented at 9 = 0, while the 
dashed line shows the phase shift for a plate with its edge at 
9 = a. 



where we have introduced the direction (i.e. angular) 
ket |f). This allows us to introduce a complete basis set 
of angular states, the eigenstates \l) of the OAM opera- 
tor C z , so that C z \l) — l\l). A Laguerre-Gaussian state 
can thus be separated into a radial and angular part, 
so that \l,p) = \pip)\l), where \pi p ) is the radial part of 
the state, which is of less importance in this paper. The 
OAM eigenstates can be transformed from integer to non- 
integer OAM states. By applying a unitary operator to 
the integer OAM basis, its completeness and orthonor- 
mality are conserved; the unitary operator that we use is 
the SPP operator, introduced in Eq. 

(f\S{\)\l) = (f|a«) = -Lexpi(l + A)0. (5) 

The new basis ||a^)| has its components enumerated 

by /, each with OAM equal to (I + X)h, where A S [0, 1) 
is a constant (not to be confused with the wavelength 
of the light). Different values of A define different bases 
of OAM Hilbert space, each basis being complete and 
orthonormal. 

As an example, passing a photon in a fundamental 
Gaussian mode |0,0) = |/?oo)|0) through a SPP with 
I = |, creates the state |/o o) l a 2/3) ' sen ding this latter 
photon through a SPP with an aligned step and I = — ^ , 
combined with a single-mode fiber, then yields a detec- 
tion probability of unity. 



V. ORIENTATION OF THE EDGE 
DISLOCATION 

Since we intend to rotate the non-integer SPPs, it is 
important to also have a description of the states that 
arise when the SPP is rotated while the coordinate sys- 
tem is fixed. To gain insight, we show in Fig.|3|the phase 
shift imposed by a SPP as a function of the azimuthal 
angle. The solid line represents the imprinted phase for 
a SPP with its edge oriented at an angle 9 — 0, while 
the dashed line corresponds to a SPP with edge orienta- 
tion 9 = a. We can now generalize the definition of the 
operator S(£) to the operator S(a, £), which includes the 
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orientation a of the edge dislocation. We find 



(f\S(aJ)\l) 



1 



/27T 



exp + £)9] 

exp [i(2n — a)£] , 
exp (— ia£) , 



< 9 < a, 
a < 9 < 2vr, 



(6) 



where a, 9 € [0, 27r). This also allows us to generalize 
the states \a^} with orientation 6 = to states with 
arbitrary orientation, |al (a)). From Eq. 10, it is im- 
mediately clear that, when neglecting the uniform phase 
shift, the complementary SPP operator S compl (a,£) is 
equal to S(a, —£) = S^(a, I). 

Since the basis ||a^'(0)) j is complete, the states after 

rotation, \aP(a)) can be written as a superposition of 

these basis states. Thus the decomposition of \a^\a)} 

into the basis ||a^(0))| depends on the angle a. To 

illustrate this, we make a projection of a non-integer state 
oriented at 9 — a, onto the same state with orientation 
6 = 0. For this, we choose a SPP with step £ = j + A, 
where j is the integer part of the step (not to be confused 
with the integer LG index I), and A € [0, 1), yielding the 
overlap amplitude 



A A (a)=<al^' ) (0)|al i+J ' ) (a)) 

=(l\&(0,j + X)S(a,j + X)\l) 
1 



2tt 



[2n — a + a exp (i2tt\)] 



(7) 



x exp [ — i(l + j + A) a] , 



where |i) is the OAM operator eigenstate with eigenvalue 
I. The overlap probability is then 

\A x (a)\ 2 = (l- ^) 2 sin 2 (A7r)+cos 2 (A^), (8) 

which depends neither on the integer part of the step 
index j nor on the OAM state I. For non-zero values of 
A, this overlap function has a quadratic dependence on 
the orientation a. This function is plotted in Fig. 0] for 
various values of A. It illustrates that, when A = 0, the 
projection does not change, as expected. For values of 
A ^ 0, the outcome of the projection is less trivial, with 
A = | providing an especially interesting case: when the 
state is rotated over a = it by rotating the SPP, the state 
is orthogonal to the non-rotated state. We will call the 
half-integer states |oi (0)) and \cSi\ir)) 'up' and 'down', 

2 2 

respectively, referring to the orientation of the edge part 
of the dislocation. In principle, this orientational label 
can be used for any non-integer OAM state, but for half- 
integer OAM states it carries an analogy to fermionic 
spin, as our 'up' and 'down' states are orthogonal just 
like up and down spin i. In polarization entanglement 
the comparable orientational labels are known as 'H' for 
horizontal and 'V for vertical polarization. To make our 
proposed OAM entanglement setup maximally equivalent 
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FIG. 4: The overlap (see Eq. (|HJ|) between a non-integer OAM 
state and the identical state rotated over a. When A = 0, the 
states are identical aside from a trivial phase shift, due to the 
vanishing edge dislocation. For half-integer OAM, i.e. A = 
the two states are generally different, leading to a parabolic 
dependence of their overlap; the two states are orthogonal 
when a = tt. 



to the polarizer setup, we focus from now on mainly on 
the case A = |. Note that, as Eq. © is independent of 
j, any SPP with £ = j + A may be used, as long as A = |. 

There is, however, a key difference between the half- 
integer OAM states on the one hand and the fermionic 
spin states or polarization states on the other hand. An- 
gularly intermediate states of the fermionic spin or of 
polarization are a superposition of the two orthogonal 
basis states 'up' and 'down', or 'H' and 'V, so that the 
overlap function varies as sin 2 a. However, from Eq. © 
we see that in the present case, the dependence of \A\ 2 
on a yields a parabola. Thus we conclude that, as the 
SPP is rotated from 'up' to 'down', the OAM state fol- 
lows a path through Hilbert space that is not confined 
to the two-dimensional subspace spanned by the 'up' and 
'down' states. Note that the OAM expectation value will 
be conserved along this path, as rotating the SPP has, of 
course, no effect on its step height, i.e. on £. 



VI. ENTANGLEMENT OF HALF-INTEGER 
ORBITAL ANGULAR MOMENTUM STATES 

When dealing with OAM entanglement, issues of con- 
servation of OAM during the spontaneous parametric 
down-conversion process arise. This conservation is dis- 
cussed in many papers 0, [H Q E E El [3 and it 
seems that OAM is indeed conserved if two conditions 
arc fulfilled, namely (i) the paraxial limit, and (ii) the 
thin-crystal limit. In practice this is usually the case and 
we will therefore work within these limits. 

The two-photon state can be described in the spatial 
polar representation as ^3 



I*) 



drT(r) \r) 1 |r) 2 , 



(9) 



where V[f) is the mode function of the pump beam. We 
restrict ourselves to a pump beam in a pure OAM mode 
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so that V(r) can be written as V(r) exp (iq8) / V^tt, where 
V(r) is the radial part of the pump mode function. Wc 
can now write the two-photon state in the half-integer 

SPP basis {|a ( x n) (0))|, thus yielding 



I*) = 



drrV(r) \r)i \r) 



E 



m.n— — oo 



d0 exp (^6»)(a ( i n) (0)|f)(a ( r ) (0)l? ; > 

2 2 

x lo^fO)) 10^(0)) 

2 2 



2tt 

where we have defined the radial ket, 

\R) = / drrV(r) |r)i |r) 2 . 



(10) 



(11) 



We now continue to calculate the coincidence fringe 
that is expected in the proposed experiment. For the 
pump beam, we shall assume a fundamental Gaussian 
beam, so that q = in Eq. (|10|l . In the signal path we 
place an analyzer consisting of a SPP with t = j + i 
where j is an integer, with its orientation set to 'up' (9 = 
0), represented by S(0,j + h), and a single-mode fiber. 
When the detector clicks, the signal state before passing 
through this analyzer is collapsed to S'(0,j + |)|0,0) = 

|poo)|cSi (0))i where |p o) is the radial part of the 

2 

fundamental Gaussian state. 

Consequently, according to Eq. 1)10(1 with I = 0, the 
idler state \ip}2 is then collapsed to 



where 



C 



|^>2 = C|a ( / ) (0)), 



((A)o|i(poo|2) l-R). 



(12) 



(13) 



Thus when analysing the idler state with S(a, j + ^)|0, 0) 
we obtain the projection 



B(a)=C(a ( P(a)\a ( P(0)) 



(14) 



Note that Eq. ((14(1 has, aside from a prefactor, exactly the 
same appearance as Eq. (JJJ) with A = \ . The coincidence 
fringe is then given by the modulus squared, 



\B{a)\ 2 = \C\ 2 (l-^-)' 



(15) 



which is proportional to Eq. (JSJ): we find a parabolic 
coincidence fringe. 



The above reasoning to obtain the coincidence fringe 
|i?(a)| 2 is valid for any choice of the signal SPP orienta- 
tion and only depends on the relative orientation a of the 
signal and idler SPPs. Thus a coincidence measurement 
on entangled OAM pairs using half-integer OAM analyz- 
ers, will bring forth a coincidence fringe that is parabolic, 
regardless of the individual settings of the analyzers. 



VII. THE CHSH VERSION OF THE BELL 
INEQUALITY 

There have been several theoretical papers that ad- 
dress the generalization of the Bell inequality [2(j 
to quantify the violation of local realism of two TV- 
dimensional particles (qu/Vits) pi EHL I2I I2I I2I : an 
example of a quAit is a spin-s particle with 2s + 1 = A. 
It has been pointed out that in this case the use of m s - 
sorting devices, such as Stern-Gerlach analyzers, does not 
offer access to higher-dimensional quantum correlations, 
presumably because the action of a Stern-Gerlach ana- 
lyzer depends only on the alignment of its quantization 
axis [2i| . Instead of larger spin values (s > ^), the use 
of spatial degrees of freedom together with Bell multi- 
ports has been advocated to gain access to the multidi- 
mensional aspects of entanglement [23|. We stress that 
all this is different from our proposed use of half-integer 
SPPs (s + i). These devices do not produce finite-A 
quAits, but imprint the infinite OAM dimensionality of 
the (oriented) edge on a transmitted light field. Rotation 
of this edge is equivalent to a partial exploration of the 
complete Hilbert space along a certain path, namely an 
iso-OAM path; due to this complexity, it is not clear how 
a generalized Bell inequality could be applied to our case. 

However, instead of using a generalized high- 
dimensional bipartite Bell inequality, it is allowed to use 
an inequality for lower-dimensional two-particle entangle- 
ment [23. Thus we choose, as in the polarization case, 
the inequality introduced by Clauser, Home, Shimony 
and Holt (CHSH) for a measurement where the coinci- 
dence probability is expected to be a function of only 
a s — a; j2(J. When relabelling a s and a; as a\ and a 2 
in no particular order), the CHSH inequality is given by 
"23 



S =E(ai,a 2 ) - E(a' 1 ,a>2) + E(ai,a' 2 ) 
+ E(a' 1 ,a' 2 ) < 2. 



(16) 



The function E is specified, for the variables x, y, as 

El 



E{x,y) 



P(x, y) + P(x ± ) y^) - P(x, y^) - P{x ± ,y) 



P(x, y) + P(x x ,y±) + P(x, J/-L) + P(x^,y) ' 

The notation x^ (and similarly for y ) is used to indicate 
an analyzer setting that analyses a state orthogonal to 
the state with setting x. Thus in our case, x^ = x + n 
and y 1 - = y + it. P(x, y) is the coincidence probability 
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function, which is equal to 



P(x,y) = \B(\y-x\)\ 2 = C 2 [l 



\y-x\ 



(18) 



As the periodicity in the present case is half that of the 
case of polarization entanglement, we use the standard 
analyzer settings for polarization entanglement 0, HH 



multiplied by a factor of two: ot\ = —\tt, a'i 



ft, 



ot2 = — tj7t, a' 2 = 0. 

Substitution yields a Bell parameter S = 3^. This is 
the key result of our paper; it indicates that in the case of 
entanglement of half-integer OAM states, the maximum 
violation of the CHSH inequality given by Eq. Ijl6(l . is 
stronger than the maximum violation that is allowed in 
polarization entanglement, namely S — 2\[2. In other 
words, quantum non-locality of the photons in the pro- 
posed setup is stronger than the maximum achievable for 
two qubits. To achieve this, only two detectors are re- 
quired and only one coincidence count rate is measured 
per analyzer setting, in contrast to the OAM quAit setup 
requiring N detectors and N 2 coincidence count rates per 
analyzer setting 0, H3| ■ 



VIII. CONCLUSIONS 

In this paper we have put forward a novel approach 
to demonstrate high-dimensional entanglement of orbital 
angular momentum states. The proposed setup uses an- 
alyzers that consist of non-integer SPPs and single-mode 
fibers, enabling detection of high-dimensional entangle- 
ment with only two detectors. 

The key idea is to use the orientation of the edge dis- 
location in the SPPs. We specialize to the case of half- 
integer I, so that the orientation of the edge as an ana- 
lyzer setting can, to a certain extent, be treated similarly 
as the axis of a polarizer in polarization entanglement. 
Instead of horizontal and vertical polarization states, we 
deal with 'up' and 'down' states, referring to the orien- 
tation of the edge dislocation. We analytically calculate 
the coincidence fringe in the entanglement setup and find 
it to be parabolic in shape, and periodic over 2tt. When 
evaluating the well-known CHSH Bell parameter, we find 
S = 3^, i.e. we predict beyond-Bell pairing of two pho- 
tons. To achieve this, we require only two detectors, as 
opposed to the standard multiport approach This 
economic exploitation of the spatial degrees of freedom 
seems to be a consequence of the singular nature of our 
half-integer SPPs, which implies, in principle, infinite di- 
mensionality. 

Experimental verification of the outlined proposal is 
under way. 
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APPENDIX: PURE EDGE DISLOCATION 

In the present paper, we have seen that, for spiral phase 
plates with half-integer step index £, i.e. imprinting a 
mixed screw-edge dislocation, the coincidence fringe in a 
twin-photon experiment is parabolic as a function of the 
relative orientation of the two radial edges, with visibility 
equal to 1. This then results in a value of the CHSH-Bell 
parameter equal to S = 3^. 

An intriguing question is whether there are other de- 
vices with optical singularities, besides the half-integer 
spiral phase plate, that, in a quantum experiment, will 
give rise to a similarly large value of the S'-parameter. 
More precisely, can one design a twin-photon experiment 
that results in a value of the S'-parameter larger than 
2\/2, using devices that are simpler to produce and han- 
dle than the spiral phase plate discussed so far? 

As we will see below, this question can be answered 
positively, and the proposed device is surprisingly simple. 
It is a non-integer step phase plate carrying a straight 
edge dislocation, as shown in Fig. Obviously, such a 
plate can be manufactured more easily than the spiral 
phase plates, whose production and characterisation is 
extensively discussed elsewhere j23, ED • 




FIG. 5: Similar to the spiral phase plate in Fig.^), the non- 
integer step phase plate has its step height h 3 chosen so that 
the phase shift due to the thick part is a non-integer multiple 
of 2n with respect to the thin part. The edge dislocation is 
chosen so that it goes through the centre of the field mode. 

Similar to the non-integer spiral phase plate, the non- 
integer step phase plate contains an orientational degree 
of freedom, a. It is possible to write a unitary operator 
for the plate's action, 



(x,y\?{a,</>)\l) 



exp (U8) 
x/2^ 



exp {iCj 
1, 



a < < a 



(19) 

where a € [0, tt) represents the orientation of the step, 4> 
the optical phase delay resulting from the step, 9 <E [0, 2tt) 
the azimuthal angle and \l) the orbital angular momen- 
tum eigenstates. For the special case </> = tt, Eq. i|19|) 
corresponds to the Hilbert transform [32l The 
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FIG. 6: The overlap (see Eq. <I21I ) between a state with an 
imprinted line dislocation, and the identical state rotated over 
a. When there is no dislocation (0 = 0), the states are identi- 
cal. For other phase differences <f), the two states are generally 
different; <j> = ^tt and (j> = n result in a continuous parabolic 
fringe. 

non-integer step phase plate obviously also transforms 
the pure orbital angular momentum states into high- 
dimensional states, due to the edge dislocation. The def- 
inition as given here can be extended for a G [— it, tt), 
seemingly complicating matters, but without any impact 
on the calculations ahead. 

Since the operator T is unitary, we can define a new, 

complete basis where T(a,4>)\l) = \b^\a)). 

We can write the states after rotation \bj (a)) in terms of 

superpositions of the non-rotated states, |&i (0)), where 
the decomposition depends on the angle a. We can il- 
lustrate this by making a projection of such a state with 
a = onto the same state a ^ 0. The overlap amplitude 
thus becomes 

A(a)=(b^(0)\b { ;\a)) 

a (20) 
=— (cos</> — 1) + 1, 

7T 

which is valid for values of a 6 [— 7r,7r). The overlap 
probability is given by 

\A{a)\ 2 = f- N ) 2 (cos0-l) 2 + 2-(cos0-l) + l. (21) 

Figure shows the intensity fringe as a function of 
a for three different values of cf>. For clarity, we show 
the functions over the range a £ [0, 2ir) instead of a E 



[— 7r, 7r), which is allowed when a is periodic over 27r. 
The most interesting cases occur when <f> = w and cf> = 
7r/2, when the fringe becomes a parabola, periodic over 
7r and 27r, respectively. The result is then, aside from 
the periodicity, identical to that for a half-integer spiral 
phase plate. 

A plate with a single edge dislocation with phase differ- 
ence tt ov it/2 can thus be used, similar to a spiral phase 
plate, to investigate entanglement. The setup would be 
identical to Fig.|21with the spiral phase plates replaced by 
the edge dislocation devices discussed here. The calcu- 
lation from two-photon state to coincidence fringe is the 
same as discussed in section lVll using the states l&i^a)), 
or \b^l 2 (a)) instead of \a^J 2 (a)). 

The CHSH-Bell inequality is violated maximally for 
these two-photon states, using the set of sixteen angles 
as used in polarisation entanglement for <j) = 7r, and the 
set as used for the half-integer spiral phase plates for 
(f> = tt/2. Both cases yield an CHSH-Bell parameter equal 
toS = 3±. 

Finally, we note that binary phase plates only slightly 
more complex than the type discussed above and pic- 
tured in Fig. can yield an even stronger violation of 
the CHSH-Bell inequality. When using a plate with the 
shape as indicated in the inset of Fig. the maximum 
CHSH-Bell parameter will be equal to S = 4, the high- 
est value at all possible. For such a plate, we predict a 
coincidence fringe with a shape as indicated in Fig. [7| 




FIG. 7: When using a plate as indicated in the inset, the 
overlap (see Eq. i'Zll ) as a function of a is as shown in the 
graph. This results in 5* = 4. 
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